Cold atom dynamics in crossed laser beam waveguides 
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We study the dynamics of neutral cold atoms in an L-shaped crossed-beam optical waveguide 
formed by two perpendicular red-detuned lasers of different intensities and a blue-detuned laser at 
the corner. Complemented with a vibrational cooling process this setting works as a one-way device 
or "atom diode". 
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I. INTRODUCTION 



Controlling the microscopic motion of atoms in gas 
phase is one of the main goals of atomic physics and atom 
optics for fundamental studies and for applications such 
as metrology, precise spectroscopy, the atom laser, or 
quantum information. Different control objectives have 
been achieved with optical and/or magnetic fields dur- 
ing the last two decades. The phase-space domain of the 
atom can be restricted by different traps, and the loca- 
tion manipulated by optical tweezers. The modulus of 
the velocity and its spread have been controlled as well 
by several stopping or cooling techniques, and its direc- 
tion by magnetic waveguides combined into atom chips 
and integrated circuits or by optical waveguides 0- 

. The implementation of complex geometries for atom 
transport is a challenging objective that may open the 
way to new interferometers and integrated quantum in- 
formation processing [ij. In particular waveguide bends 
are basic elements that have been investigated experi- 
mentally and theoretically [6l-[l0|. 

Aside from modulus and direction, the control of the 
remaining element of the atom velocity as a vector, 
its sense or orientation (say to the right or left for a 
given direction), has been undertaken much more re- 
cently with theoretical proposals and experimental pro- 
totypes of atomic one-way barriers or "atom diodes" [lll - 
[23j . They are analogous to a semipermeable membrane 
or a valve, which let the atoms cross it one way (for- 
ward) and block their passage in the other one (back- 
ward). A conceptual precedent is the automated demon 
conceived long ago by Maxwell to achieve a differential 
of pressure between two parts of a vessel and demon- 
strate the statistical character of entropy IJ] . Maxwell 



only specified the demon's action, not its inner workings, 
whereas, more than one century later we are beginning to 
design and realize such devices. Applications that have 
motivated so far this research are the possibility to cool 
species without cyclic transitions (igI . [22j . or the con- 
struction of trapdoors and flow control in atomic chips 
and circuits [ll|. The existing methods are essentially 
one-dimensional (the controlled sense corresponds to a 



longitudinal or a radial velocity). A basic scheme con- 
sists on setting a barrier in one atomic level, e.g. the 
ground one. On one side of the barrier, say the left, the 
atom is excited adiabatically so as to avoid the ground 
state barrier. Adiabaticity is useful to make the transfer 
efficient and velocity independent in a broad range and 
to avoid the passage from the upper to lower level for 
atoms that approach the diode from the left [13 ■ On 
the other side of the barrier the excited state is forced to 
decay so that an atom coming from the right is reflected 
by the barrier. One possible variant is to substitute the 
adiabatic step by optical pumping. It is also velocity 
independent in a broad range and precludes forced deex- 
citation on the wrong side. An irreversible step, which 
ideally may be reduced to the emission of one-photon, 
is essential to break time-reversal invariance, a necessary 
condition to realize a true one-way barrier. Any atom 
diode has of course certain limitations with respect to 
velocity working range, efficiency, width of the structure, 
or species and states that can be treated. For example, in 
a two-laser, optical prototype for Rubidium [l^, Hj], the 
barrier produced some undesired heating because the in- 
ternal hyperflne structure used does not allow for a suffi- 
ciently large detuning; the use of magnetic sublevels may 
avoid this effect but at the price of loosing many atoms 
because of the branching ratios during optical pumping 
[T6| . While these limitations may or not be relevant de- 
pending on the intended application, it is desirable to 
investigate other mechanisms, surely subjected to differ- 
ent constraints. 

In this paper we investigate two aspects of cold atom 
guiding and their possible combination into a single de- 
vice: bends in L-shapcd asymmetric guides and one- 
way motion. Most previous studies of bent waveguides 
have focused on magnetic implementations and sym- 
metrical arms @, l8l-[l0|: AT-shaped optical waveguides 
have been investigated as beam splitters for interferom- 
etry 0, US t25j. An experimental realization of an X- 
shaped asymmetrical beam splitter, with different po- 
tential depths in the two guides, has been also carried 
out [1^. We shall study here an optical, asymmetric re- 
alization of an L-shaped bend and determine the trans- 
fer between longitudinal, gap, and transverse energies. 
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an angle clockwise with respect to the y-axis, more on 
this below. The corresponding potential is 



100 



FIG. 1: (Color online) Scheme of the potential created by a 
blue detuned and two red detuned lasers. U' = 0.474 /xK. 



In addition, when combined with vibrational cooling the 
L-shaped guide provides a 2D one-way mechanism for 
one-way motion. 



II. SIMPLE OPTICAL WAVEGUIDE MODEL 

The proposed device consists of three horizontal Gaus- 
sian laser beams (see Fig. [T]). Two of them are detuned 
to the red with respect to a transition between the atomic 
levels g, e and play the role of waveguides for the ground 
state atoms along the x and y axes. We shall perform 
2D simulations corresponding to a tight confinement in 
z (vertical) direction, ignored hereafter, by an optical 
lattice. We have simplified the corresponding potentials 
neglecting the dependence with the longitudinal coordi- 
nate, 

U{x,y) = -Ue-^y'/""^ > 0, upper valley, (1) 
U'{x,y) = -t/'e^^"^'/"'" t/' > t/, lower valley. (2) 

(w and w' are the waists.) This is reasonable within 
the Rayleigh length. It is as well a simplified treat- 
ment for combined magneto-optic waveguides Q in which 
the longitudinal potential dependence is essentially sup- 
pressed by cancellation between a repulsive magnetic po- 
tential and an attractive optical potential. Note that the 
assumed asymmetry in intensities creates "upper" and 
"lower" valleys in the potential energy surface.^ 

A third laser, detuned to the blue, forms a barrier to 
redirect the atoms from the upper to the lower valleys 
blocking the passage to the red detuned arms along the 
positive-x and positive-y semiaxes. This laser is displaced 
slightly away from the coordinate origin and it is rotated 



"Upper" and "lower" refer to the energy, not to a relative spatial 
height. Quantities such as energies, velocities or momenta associ- 
ated with the upper/lower valley will be umprimmed/primmcd. 



lJ^{^x,y) = c/(je~^[^'^~^'')cose+(a-ao)sine]^/tog 



(3) 



We shall study the atom dynamics with quantum ap- 
proaches (wavepackets and stationary methods), and 
with classical trajectories too since they provide a rather 
accurate description -in particular when an average over 
the transverse phase is performed- in a much shorter 
computation time than the quantum simulations. For a 
given incident longitudinal energy and vibrational state 
we do not perform "Ehrenfest" (one trajectory) classical 
simulations but ensemble averages over all possible 
phases of transverse motion to avoid the sensitivity of 
classical trajectories with respect to the phase and bet- 
ter mimic the quantum results. The details are given in 
Appendix El 

We assume that there is no significant interference 
among the three beams so their potentials simply add 
up. This may be achieved e.g. by orthogonal polariza- 
tions of the red detuned lasers and/or different detunings 
that cause a fast time-dependent interference that aver- 
ages out in the scale of the atomic motion [27]. 

In wavepacket computations, see Appendix |B] for nu- 
merical details, we assume that the wave function of the 
initial state factorizes into longitudinal and transversal 
functions, 



^'(a;,y,0) = V(a;,0)(»4'(2;,0). 



(4) 



For atoms incident in the upper channel the initial trans- 
verse wave function $(y,0) will be the ground state of 
the Gaussian potential, Eq. ([T]), which is calculated nu- 
merically by diagonalization of the Hamiltonian. In the 
longitudinal direction we choose a minimal uncertainty- 
product Gaussian, 



V'(a;,0) - 



1 



(2^^2)1/4 



(5) 



where is the width of the wavepacket and x^ and pi 
the initial position and average longitudinal momentum 
respectively. For atoms incident from the lower channel 
a corresponding approach is used interchanging x and y. 



III. FORWARD MOTION: PASSAGE FROM 
THE UPPER TO THE LOWER VALLEY 

We shall discuss first the main factors that determine 
the passage of atoms from the upper valley, incident in 
the ground vibrational state, to the lower valley. All 
calculations are done for the mass of ^'^ Kb atoms. 

The harrier position. If the barrier is too far from the 
crossing point of the waveguides, a well is formed due 
to the addition of the upper and lower valley potentials, 
Eq. dl]) and Eq. ([2]), see Fig. Wp- This well allows for 
long lived chaotic (classical) trajectories and favours en- 
ergy transfer among the degrees of freedom as well as 
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FIG. 2: (Color online) Potential contour maps and representa- 
tive trajectories for different wall positions, (a) xo = yo ~ 0; 
(b) Xo = yo = 5/im; (c) xo — yo = 3^m. Laser param- 
eters: U = 0.158 fiK; U' = 0.474 fiK; Ut = 0.948 ^K; 
w = w' = wt ^ 5.7 ^m and 6* = 45°. The 3D plot of Fig. 1 
corresponds to case (c). 
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FIG. 3: (Color online) Classical and quantum probability to 
find the atoms in the lower valley for different rotation angles 
of the blue detuned laser. Wavepacket parameters: = 2 
fim; Xi = —15.88 fim; Vi = pi/m = 0.41 cm/s. (For the 
classical calculation we average the fixed-energy probabilities 
with the momentum distribution of the longitudinal Gaus- 
sian.) Laser parameters: U — 0.158 fiK; U' — 0.474 ^K; 
Ub = 0.948 /iK; w = w' = = 5.7 pim and xo = = 3 /im. 



reflection back into the upper valley. Displacing the blue 
detuned laser nearer to the origin the well is filled and 
the chaotic behavior and reflection are avoided. The wall 
should not be too close to the crossing though, as it would 
obstruct the upper valley and thus the atom passage, as 
in Fig. [2^. Between the two extremes there is a range 
of distances for which the well is suppressed without ob- 
structing the upper valley, see Fig. [2t. Representative 
classical trajectories for the three cases are depicted in 
Fig.H 

Barrier angle. Remarkably, the probability to pass 
from the upper to the lower valley shows a stable full- 
transmission plateau for a broad range of angles 0. This 
is shown in Fig. [3] for wavepacket and classical trajec- 
tory calculations. The optimal choice of angle depends 
on its effect on energy transfer among longitudinal and 
transverse degrees of freedom, as discussed next. 

Vibrational excitation. If the atom passes to the lower 
valley, the asymmetric potential configuration favors its 



FIG. 4: (Color online) Average final vibrational energy versus 
initial kinetic energy for different crossed-beam setups com- 
puted with classical trajectories. The dashed line is a quan- 
tum calculation (for monochromatic incident energy) showing 
good agreement with the classical one. Solid line: K — {V). 
Laser parameters: w = w' = 1.03 /im, wt = 0.77 /im, 
a;o = yo = 0.41 fim, U = 0.4 mK and Ut ~ 5.2 mK (rombs). 
Rest of cases: w = w' — wt = 5.7 fim; xo ~ yo = S ^m; 
U = 0.158 fjK, Ub = 0.948 fiK (squares, triangles and stars), 
and Ub = 1.9 fjK (circles). The lines end when the transmis- 
sion probability is no longer one. This full-transmission range 
increases with the angle and the depth of the lower valley. 



vibrational excitation (or "transverse heating"). For a 
transition from the ground state of the upper valley (rt = 
0) to the n' vibrational state of the lower valley (0 — >■ 
?i' for short) conservation of energy, measured from the 
bottom of the lower valley, takes the form 



E = K + Vo 



(6) 



where A = — is the gap between valleys, K, K' 
are the upper and lower kinetic energies, and Vq, Vn' the 
corresponding vibrational energies (measured from the 
bottom of each valley). Fig. 2] shows the (lower valley) 
average vibrational energy {V') versus the incident K 
for several cases. Even for K ^ Q the process is highly 
non-adiabatic (a simple ID adiabatic treatment as in [3] 
is therefore not valid here), and a significant fraction of 
the potential energy gap is converted into vibrational en- 
ergy. As K increases, the trajectories penetrate more 
on the reflecting blue wall so that the outgoing trajecto- 
ries are further away from the bottom of the lower valley 
and vibrational excitation increases. The average {V') 
is essentially linear in at variance with a quadratic 
dependence found for circular bends @. At the bottom 
of the lower valley the kinetic energy of a classical trajec- 
tory equals the total energy iiT -I- A -I- Vq (measured from 
the bottom of the lower valley). It may be split into x 
and y components taking into account the angle a of the 
velocity with the y axis. The cc-component is the final 
vibrational energy and it takes the form 

V' = Esm^{a) = {l^ + Vo) ain^ [a) + K sin^ {a) , (7) 

but a is roughly constant for a given set of potential 
parameters because of the relative flatness of the impact 
region at the waveguide corner. This region results from 
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FIG. 5: (Color online) The quantum transmission probabil- 
ities obtained from Eq. (|B9|) for different vibrational lev- 
els of the lower guide versus the initial velocity. Laser pa- 
rameters: U = 0.158 ^K; U' = 0.474 ^K; Ut = 0.948 fiK; 
w = w' = Wb ~ 5.7 /im; xq = ya — 2i /im and 6 — 45°. 



the combination of the dominant lower valley and barrier 
potentials. 

Most lines in Fig. 4 are for classical-trajectory compu- 
tations but we have also checked the good agreement with 
a fully quantum calculation in one case. To do so we have 
extracted the quantum, stationary (fixed energy) state- 
to-state transition probabilities — > n' , {qn' /p)\Ton' 
from wavepacket calculations as explained in the Ap- 
pendix [B] Here g„' is the longitudinal momentum in the 
lower valley for the vibrational state n' and Tqu' (p) the 
transmission amplitude for incident longitudinal momen- 
tum p = \j2mK. In Fig. [S] we show the dependence of 
the quantum transmission probabilities versus the initial 
velocity. Note again, now in more detail, the increase of 
vibrational excitation with p. For sufficiently large en- 
ergy this leads to escape from the trap. 

Fig. [S^ shows the total transmission probability Pt = 
S^=o ^l^0n'(p)P versus initial velocity for 9 = 45°, 
N' being the maximal vibrational number in the lower 
valley. Note the good agreement between the quantum 
and classical calculations. Pt is very stable, and only de- 
cays from one due to escape from the waveguide caused 
by the increasing transverse heating. In principle the en- 
ergy threshold for escape is, from conservation of energy, 
if + Vo + A > U' (solid vertical line on Fig. but 
the effective threshold occurs at higher energies, when 
ly') K, U' , since the available initial total energy is trans- 
ferred only partly into vibrational energy (V); a sign of 
the escape is the coincidence of decay of Pt with the pop- 
ulation of the highest vibrational level {N' = 40 for the 
chosen parameters). The velocity range for full forward 
transmission may be increased at will, according to Eq. 
([7]), by increasing the gap A, an example is shown below. 

A second effect that may spoil the forward passage is 
the possibility to overcome the barrier when K + Vq > 
Ub (we neglect here the lower valley potential). This 
threshold is higher than the former, and is marked by a 




E/U' 



FIG. 6; (Color online) (a) Classical and quantum probabil- 
ities to find the atoms in the lower valley when they start 
in the n = vibrational level of the upper valley, versus the 
initial longitudinal velocity. The vertical lines are the energy 
thresholds to escape from the guides (solid) and overcome the 
barrier (dotted), (b) Classical probability to find the atoms 
reflected in the lower valley when they start in the n' = 
vibrational state of the lower valley versus the initial veloc- 
ity {Vn' = <ln' /fn). Vertical lines mark the energy thresholds 
to pass to the upper valley (dashed), escape from the guides 
(solid), and overcome the barrier (dotted), (c) Combination 
of (a) and (b). Blue circles: Transmission probability P = Pt 
versus the total energy (measured from the bottom of the 
lower valley), the atoms start in the upper valley. Red trian- 
gles: P = 1 — Pr versus the total energy of the sample, the 
atoms start in the n' = vibrational level of the lower val- 
ley. Dashed arrow: It indicates the total energy that remains 
for backward motion after the (perfect) vibrational cooling 
process. Laser parameters: U = 0.158 /xK; U' = 0.474 /xK; 
Ub = 0.948 /xK; w = w' = wt = 5.7 nm; xo = yo ~ 3 /^m and 

e = 45°. 



vertical dotted line in Figs. and [5^. These two effects 
are illustrated with representative trajectories in Fig. [71 
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IV. OBSTRUCTED PASSAGE FROM THE 
LOWER TO THE UPPER VALLEY 
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FIG. 7: (Color online) Classical trajectories that illustrate 
two guiding breakdown mechanisms: escape by vibrational 
excitation (dashed line v = 1.5 cm/s), and surmounting the 
barrier (solid line, v = 1.8 cm/s). Laser parameters: U — 
0.158 fiK; U' = 0.474 /iK; Ut = 0.948 ^K; w = w' = wt, ^ 5.7 
fj,m; xo ~ yo ~ 3 /ini and 9 — 45° . 



The potential asymmetry causes an asymmetry in the 
dynamics since, in general, for the same total energy, the 
probabilities Pq' and Po-yn'(-E') are quite different. 
This is compatible with time reversal invariance, which 
imphes only the equality for probabilities of a process and 
the time reversed one. P„^„'(£') = Pn>^n{E) holds as 
long as no irreversible step takes place (that case will be 
considered in the following section). The nature of the 
stated asymmetry can be understood from the potential 
contour in Fig. [Jt, or the 3D plot in Fig. [T] Even 
when the passage n' = — >■ n = (0' — ?> for short) 
is energetically allowed, a vibrationally unexcitcd atom 
does not find easily the lateral gate to the upper valley so, 
for a range of energies above the energy threshold (dashed 
vertical line in Figs. [Sb and|H}D), the atom is still reflected 
into the lower valley. This may be seen in Figs, and 
[5Jd, where the reflection probability is shown for states 
beginning in the fundamental vibrational state n' = of 
the lower valley for 9 — 45° and 9 — 60°. The dynamical 
reflection is enhanced by increasing the angle 9 so that 
the backward collision is more head-on, but increasing 
it too much may obstruct the passage in the forward 
direction at low velocities. Above the energies with full 
reflection, the atom with backward incidence may escape 
from the guides, when Vn' + K' > U' , or surmount the 
potential barrier when Vn' + K' > Ub (in this inequality 
we neglect the small effect of the upper valley potential) . 
The corresponding energy thresholds for these processes 
are marked by solid and dotted lines in Figs. and[HjD 
but, as for forward motion, the effective thresholds occur 
at higher velocities. 
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FIG. 8: (Color online) (a,b,c): Same as in Fig. [6]for a different 
laser configuration. Laser parameters: U — 0.4 mK; U' — 2.8 
mK; Ub — 5.2 mK; w = w' = 1.03 ^m; wt = 0.77 ^m; 
xo ~ yo = 0.41 /im and 9 = 60°. 



V. DIODE EFFECT 

The stable plateaus for full transmission and reflection 
and the asymmetry for forward and backward motion 
from the ground transverse states are prerequisites for 
a diode but not enough. A diodic or "one-way" barrier 
effect is achieved by complementing these features with 
vibrational cooling in a region of the lower valley. Several 
cooling mechanisms have been demonstrated or proposed 
for neutral atoms in tight traps: Tuchendler et al. [1^ 
have cooled single ®''Rb atoms in the tight- confining di- 
rections of a strongly focused dipole trap with optical 
molasses; Sideband cooling has been demonstrated for 
alkali-earth atoms using a "magic weavelenth" light- 
shift compensating technique [13], and for Cs atom s by 
means of 2-photon Raman transitions in ID 'sif, 2D [32| . 
and 3D [33] far-detuned optical lattices; rf-induced Sisy- 
phus cooling has been also realized for ®''Rb [33 ]. 

We shall not model in detail any of these methods here 
but simply assume that vibrational cooling is performed 
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on the atoms that have been heated transversely in the 
forward passage and analyze the consequences for back- 
ward motion. In the ideal case of cooling down to the 
ground state, n' — 0', keeping the same kinetic energy 
K' , the backward passage to the upper valley is energeti- 
cally forbidden if K' + Vq' <Vq + A, or using Eq. dH) and 
replacing K' by (V), (V) > K (Vb' is neglected). In fact 
it will not occur even at higher kinetic energies because of 
the 2D reflection effect described in the previous section. 
To determine if backward reflection is possible for a given 
incident if, we need the forward transmission and back- 
ward reflection velocity intervals of Figs. [n^,b or |5^,b 
(this has been combined in Figs. [6]: and [S};, where the 
reflection information is represented by 1 — -Pr), and the 
dependence {V'){K). After (perfect) vibrational cooling 
the backwards energy is 

K' + Vo, ^K + A + V^-{V')+V^,, (8) 

We may now check the value of Pr{0') for this energy 
to see if the atoms are reflected back into the lower val- 
ley. We have done this for the edge points of the total- 
transmission interval and the result is represented by the 
arrows in Figs. |BJ: and IHt- Note that for the 45° case in 
Fig. IB}:, the high velocity edge of full forward transmis- 
sion does not correspond to backward reflection. This 
may be remedied by increasing the range of full reflec- 
tion with a larger angle. For a lower valley as deep 
as the trap in [28[ and 9 = 60°, see Fig. [B}:, a broad 
stable operating range for diodic behaviour is achieved, 
where the full range of forward passage corresponds, af- 
ter transverse cooling, to full reflection in the backward 
direction. 



VI. CONCLUSIONS 
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Appendix A: Classical dynamics 

Classical trajectories are a useful tool to explore the ef- 
fect of varying parameters faster than the quantum com- 
putation. They also provide physical insight. We solve 
Newton's equations 

dU{x,y) .. dU{x,y) 

mx = , my = , (Al) 

ox ay 

where U{x,y) — U + tj' + Ub are given in Eqs. (HI [21 
[3]), transformed into a system of four equations with a 
fourth-order Runge-Kutta method. 

To mimic the scattering at fixed longitudinal and vi- 
brational energies we consider an ensemble average set as 
follows: we take first a classical reference particle mov- 
ing periodically in the transversal direction of the upper 
valley with the same transverse energy as the quantum 
state. To run an even number N of trajectories the pe- 
riod of this reference particle is divided into N equal time 
segments and the N values of ti set the initial 

transverse conditions for the trajectories of the ensem- 
ble. For the longitudinal motion we simply impart to the 
trajectories the longitudinal momentum p. 



Guided atom lasers in the ground state of the trans- 
verse confinement have been recently realized 0-111 and 
more complicated settings are being considered, in partic- 
ular with crossed beams, following similar developments 
in magnetic waveguides that may pave the way to new 
interferometers, atom integrated circuits and analogs of 
electronic devices [H, [sl] . 

In this work we have explored a realization of straight 
angle bends in asymmetrical optical waveguides for cold 
atoms, with two red and one blue detuned lasers, as well 
as the possibility to use the transverse heating caused 
by this geometry, combined with vibrational cooling, to 
implement a diodic (one-way) device. Indeed the trans- 
mission and reflection probabilities of the proposed struc- 
ture offer the stability with respect to incident velocity 
required for an efficient diode. The different elements 
of the proposed device have been already implemented 
separately, and the remaining technical challenge is their 
combination into a single device. 



Appendix B: Split-Operator Method (SOM) 

Given the time dependent Schrodinger equation with 
Hamiltonian 

H=|l+U=4T^ + t>(^,n, (Bl) 
Am Zm 

the Split Operator Method (SOM) approximates the evo- 
lution operator as 

The resulting integrals are easily solved using the Fast 
Fourier Transform (EFT) technique [37| . 
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1. Discretization and experimental setting 

The validity of the discretization approximation re- 
quires [13, HI 



Q 



1 



Ax 

Tlx Lx ^TtQol 

hn 

da 

h h 

dt < 



Pjtl.x 7 ) 

ax 



T ' TI 

J- max ^ rt 



(B3) 
(B4) 
(B5) 



where Q^: is a quahty factor which takes into account 
the number of the lattice points that represent the wave 
function in coordinate and momentum representations, 
Lx and rix are the lattice length, and the number of di- 
visions, and dx and dt are the space and time steps. Ax 
and pm are the minimal spatial dispersion (usually the 
one at t = 0) and the maximum momentum value. Fi- 
nally Umax is the maximum potential energy and Kmax 
is the maximum kinetic energy during the simulation, 

Kmax,x < 2m " ^'^'^ y-dircction we have similar con- 
ditions. In all calculations we set Q > 15. Other param- 
eters are n-x ^ Uy = 4096. 



2. Stationary transmission amplitudes from 
wavepacket computations 

We write the transmitted wavepacket state as 

/oo 
- OO 

n- 



where p is the incident, longitudinal momentum of the 
atoms in the upper valley, q„/ — [p^-f2m(Vo-|- A— y„/)]^/^ 
is the longitudinal momentum for the vibrational state 
n', Vn'{x) — {x\vn') IS the amplitude of a lower valley 
vibrational state and ipip) = {p\''Pi^j 0)) is the initial mo- 
mentum distribution of the wave function given by Eq. 
([5]). Finally the transmitted wave function is projected 
onto one particular eigenstate Vn' (x) , 



{vn'\'i/T{x,y,t)) = / dp Ton'{p)(j){p) 



.(B7) 



Defining the inverse Fourier transform as 



K'|*T(x,2/,t))e*^*/\ (B8) 



and integrating Eq. (jB7|) with respect to time from — oo 
to oo (In practice t — Q plays the role of t = — oo, whereas 
i = oo is approximated by the time when the tails of 
the transmitted wave function are not affected by the 
barrier.), we obtain 



7o„' {p) 



p e 



iy^p'^+2m{/\+Vo-V,^,)/h 



m 



2m 



(B9) 

Its modulus squared times qn' jp gives the transmittance 
(transmission probability) from the ground state of the 
upper channel to the n'th-vibrational level of the lower 
guide. 
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